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Abstract 

Precision measurement of the leptonic CP violating phase 5 will suffer from the, then surviving, 
large uncertainty of sin 2 #23 of 10 — 20% in the experimentally interesting region near maximal 
mixing of #23- We advocate a new method for determination of both #23 and 5 at the same time using 
only the v e and v e appearance channels, and show that sin 2 #23 can be determined automatically 
with much higher accuracy, approximately a factor of six, than sin 5. In this method, we identify a 
new degeneracy for the simultaneous determination of #23 and 5, the #23 intrinsic degeneracy, which 
must be resolved in order to achieve precision measurement of these two parameters. Spectral 
information around the vacuum oscillation maxima is shown to be the best way to resolve this 
degeneracy. 
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I. INTRODUCTION 



In neutrino oscillation physics, the traditional way of determining the leptonic CP vio- 
lating phase, 5, is to measure 6 together with $ 13 by using u e and u e appearance channels 
in muon neutrino superbeam experiments, (or similarly by and appearance in either 
neutrino factory or beta beam experiments). In posing the problem in this way, $23 is as- 
sumed to be determined to high accuracy by a disappearance measurement up to the 
octant degeneracy. But, now we need a new orientation to choose the right approach for the 
determination of 5, because: 

• $13 is already determined with high accuracy, and its precision will become even greater 
by the time 5 is measured 

• and $23 will be determined with less percentage accuracy than $ 13 in the experimentally 
preferred region, see (1,0], between 40° and 50°. Since the disappearance measure- 
ments will have difficulty in determining S23 = sm2 $23 with an accuracy better than 



How seriously does the uncertainty in affect the sensitivity to 51 We show in Ap- 
pendix |X] that the uncertainty of 5 is approximately proportional to 



where P and P are the v e and v e appearance probabilities with uncertainties AP and AP, 
respectively, and As^ is the uncertainty on s| 3 , and 77 is generically a number of order unity. 
Thus, when the measurement enters into a precision era in which the percentage uncertainty 
on P and P would be smaller than, say 5%, the uncertainty of 5 would be dominated by 
the S23 uncertainty 1 . 

In this paper, we present a new strategy for precision measurement of 5 and $23 to 
overcome this problem. Instead of using disappearance channel for determination of $23, 
we rely on v e and v e appearance measurement to determine precisely $23 and 5 at the same 
time. With this method the uncertainties expected for $23 and 5 are related with each other 
as follows 



near the first vacuum oscillation maximum, as shown in Sec. IIII Bl Therefore, once we 
enter into the era in which sin 5 can be measured with reasonably high accuracy, S23 can 
be determined automatically with significantly higher accuracy than sin 5 which can be 
competitive with the disappearance measurement in the region 40° < $23 < 50°. Of course, 
the disappearance measurement will provide supplemental information, and also serves as a 
consistency check on z/SM. 

This new setting, of simultaneous measurement of $23 and 5, brings us to the problem of a 
new degeneracy involving these two parameters which we call the "$23 intrinsic degeneracy" 
because of it's similarity to the $ 13 intrinsic degeneracy [§]. Notice that this degeneracy is 



Note that this is a simplified discussion for the sake of clarity. A more detailed and accurate description 
can be found in Appendix A. 
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entirely different from the #23 octant degeneracy discussed in [lGj. We show that resolving 
this #23 intrinsic degeneracy is crucial to achieve precision measurements of 5 and 623, and 
that the disappearance measurement cannot help in region sin 2 2623 ^ 0.97. We argue that 
spectrum information around the vacuum oscillation maximum peaks is the most powerful 
way to resolve the degeneracy. It may be possible that such spectrum analyses could be 
carried out with a narrow-band neutrino beam, but if it is not powerful enough, use of a 
wide-band beam could be required. 

If the mass hierarchy is unknown the #23 intrinsic degeneracy would be a part of a larger 
degeneracy similar to the conventional #13 — 5 eightfold parameter degeneracy . Here, 

we focus on the "intrinsic" part by simply assuming that the mass hierarchy will already be 
determined at the time when the precision measurement we discuss will be realized. The 
generalization to the case including the unknown mass hierarchy can be performed in a 
straightforward manner. 

This paper is organized as follows: in Sec. II we introduce the 823 intrinsic degeneracy, 
the stimultaneous measurements of 5 and 623 is discussed in detail in Sec. Ill & IV. In Sec. 
V & VI a simple toy analysis is given for a variety of possible future precision experiments. 
The conclusions are presented in Sec. VII. The three Appendices contain some mathematical 
derivations too detailed for the main body of the paper. 



II. INTRODUCTION TO THE 6 23 INTRINSIC DEGENERACY 

We start by giving pedagogical reminder of the relevant issues, some of which have been 
discussed in the past, but not all of them. Suppose one could measure sin 2 2# 2 3 to an 
uncertainty of 0.02 at, say, 95% C.L. Then, to what accuracy can one determine sin 2 6*23 = 
S23? The answer to this question is given in the left panel in Fig. [TJ In the region of 623 
far away from maximal, the loss of sensitivity to s| 3 is modest, apart from the octant issue, 
whereas for sin 2 2# 2 3 ^ 0.98 the accuracy of s| 3 suddenly jumps to 10-20% level. This was 
observed in Fig. 3 of ref. [s| where they plotted the s^-error as a function of sin 2 2#23- 

In the left panel in Fig. [TJ two vertical scales are provided for convenience of the readers. 
Observe that the nonlinearity of scales in mapping from S23 to sin 2 2# 2 3 is large, and it is 
most significant in region where 

sin 2 2fl 23 > 0.97 or 0.41 < sin 2 fl 23 < 0.59 or 40° < fl 23 < 50°. (3) 

This is the region with significant overlap with the experimentally preferred one 0, 0] , and 
it is the region we are primarily interested in for the purposes of this paper. 

In the right panel of Fig. [TJ the bi-probability plot [11], gives the allowed values of 
(P(^ u e ), P(z/ M P e )) for three set of values of sin 2 29 23 =0.9516, 0.9676 and 0.9996. 
Here, we have switched off the matter effect and restricted ourselves to the Normal Hierarchy. 
The figure shows that with an exact measurement of the v e and v e appearance probabilities 
one can determine s\ z and 5 at the same time apart from the twofold degeneracy, which 
we call the 623 intrinsic degeneracy. Notice that this degeneracy is completely different in 
nature from the 623 octant degeneracy and exists for all values of 6*23, even at 7r/4. 

To illuminate the nature of the new degeneracy on 823 and 5, we place in the left panel 
of Fig. [T] three pairs of cross marks, distinguished also by color, corresponding to three sets 
of ellipse shown in the bi-probability plot, the right panel. In the two cases far off maximal 
mixing 6*23 depicted by the blue and the red crosses, it is obvious that the degeneracy has 
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FIG. 1: Left panel: Allowed region in sin 2 2023 vs - sin 2 #23 space assuming 2% uncertainty of 
sin 2 2#23 by a disappearance measurement. Right panel: The bi-probability plot in Piv^ — > v e ) 
vs. P{pii — > v e ) space with three set of values of Piy^ — > Ve) an d P{pn — >■ v e ) which correspond to 
sin 2 2^23 approximately 0.96 (blue and red), 1.0 (black). The blue and red pair of ellipses are in 
different octants. The precise value of sin 2 #23 of each ellipse is as indicated in the figure, and at 
the crossing points indicated by x the value of 5 is 60° (120°) for larger (smaller) sin 2 #23 ellipses, 
for the normal hierarchy. The colored crosses on the left panel correspond to the ellipse on the 
right panel. 

nothing to do with the #23 octant degeneracy. In the case near the maximal 823 drawn by the 
black crosses in Fig. [H the two solutions live in different octants, and we have a merging of 
the #23 intrinsic degeneracy and the octant degeneracy for 823 at very close to the maximal 
mixing. 2 

We note that the size of the ellipses does not change much with different values of 823 
because they scale as oc sin 2#23 which changes very little near 823 = vr/4. If the value of 
813 had been changed instead of 823, not only does the location of the bi-probability ellipse 
change but also its size changes with increasing or decreasing sin #13. Note also that in 

vacuum with arbitrary values of A = A ™g L , and in matter at vacuum oscillation maxima 
(VOM, A = (| + nj 7r), the 6*23 intrinsic degeneracy is a 5, n — 5 degeneracy. It is because the 
difference between the degenerate solutions of sin 5 vanishes in vacuum (see (j!4p ). and only 
sin 5 is determined at VOM. This is reminiscent of the 813 intrinsic degeneracy first discussed 
in [9| but there are important differences. Understanding the nature of this degeneracy will 
be elaborated in the following three sections, IHIt HV] and IVl 

From the bi-probability plot given in Fig. [H one can infer an interesting feature of the 5 
dependence of the error of 5. The uncertainty of 823, which allows the ellipses moving up 
and down, results in larger errors for 5 determination at the both ends of the ellipses, i.e., 



2 If the precise value of #13 is known, adding only disappearance measurement produces the #23 octant 
degeneracy. Adding only v and v appearance measurement produces the #23 intrinsic degeneracy, whose 
solutions may exist across the octant boundary for #23 very close to maximal mixing (black line and 
crosses in Fig. [IJ. With a finite error of #13, the two degeneracies fuse in this region. 
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around 6 ~ ±f , than at 5 ~ 0, or ir. Therefore, the large uncertainty of s| 3 described in (J3]) 
can produce broad peaks at 5 ~ ±| in the 5 dependence of <5 error. In fact, this feature is 
seen in some of the figures (e.g., in Fig. 9) given in (l2| . 



III. CORRELATED MEASUREMENTS OF 
ANCE CHANNELS 



7 22> 



AND 5 BY v P AND v P APPEAR- 



We demonstrate our new strategy by showing that simultaneous determination of #23 
and 5 is the way to proceed. To illuminate the point we consider a counting experiment to 
measure the v e and v e appearance probabilities at a given energy in matter; 
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A® is of order unity in most regions of relevant experimental parameters, but A e and A e 
are suppressed by smallness of S13. The A's are defined such that the atmospheric and the 
solar oscillation probabilities at the maximal mixing angle #23 = f , P® and P & respectively, 



are given by P e = A® and P Q = e A Q . For simplicity we will ignore the 2c 23 e A & terms, in 
Eq. (jlj), in the rest of this paper as they are significantly smaller than the other terms. 3 

It is instructive to quote an estimation of the ratio of A® to eA Q in near vacuum environ- 
ment 
A ~ 



<C 1, assuming the experiment is performed near the first oscillation maximum, 
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As we will show later, the size of this ratio is significant for the relative size of the uncer- 
tainties for the measurement of sin 2 9 23 and 5. 



3 To include the effect of 2e 2 c l , 3i A 2 V) terms, rewrite c 2 R = 1 — s| s and then the 2s 2 , R e 2 A% ] can be included in 
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the A\ terms and the 2e 2 A%, can be moved to the left hand side. 
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A. Measuring #23 and 5 at the first vacuum oscillation maximum 



To show the point let us consider measurement at the first vacuum oscillation maximum 
(1st VOM), A = |. Then, cos (5 + A) = + sin 5. One can eliminate 5 by combining the two 
equations in (ffl, then 
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(8) 



If instead we eliminate the first terms of flU) in favor of sin 5 we obtain, 



sin<5 = 



2 sin 2# 23 (A ffi + A e ) eA Q 
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(9) 



cos 5q can be obtained as cos 5q = + y 1 — sin 2 5q. Notice that the sign ambiguity 5 <H> n — 5 
cannot be resolved by measurement here. Apart from the different linear combinations of P 
and P, these two equations differ by the overall factors of 1/ y A^A^ and 1 / eA Q in s\ 3 and 
sin 5, respectively 

We use the matter perturbation theory to obtain a simpler expression which allows in- 
tuitive understanding. We assume .% <C 1, which is an excellent approximation for T2K, 
then we obtain, to first order in 
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(10) 



where J r = C12S12C23S23S13C13. Therefore, at VOM and in near vacuum environment s| 3 is 
determined by sum of the probabilities P + P, while sin S is governed by the difference P — P. 
P + P and P — P switch their positions to induce sub-leading corrections due to matter 
effect into sl 3 and sin 5, as shown in ( ITUl) . 



B. Error of is smaller than that of sin S 

Let us estimate the expected errors in measurement of s| 3 and sin 5, assuming that the 
errors for P and P are independent with each other. We also assume settings in which the 
matter effect remains subdominant and use notations A™ c for A^ = in vacuum. Using 
([8]) and ([9]), the errors, A(s| 3 ) and A(sin<5), are given by 



A (40 - 77lkW(^) 2 + UP)> 



1 



Notice that the two errors are related with each other at the 1st VOM by 

A (4>) « (|^|) A ( sin ^) - ^ A ( sinJ ) ( 12 ) 
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Therefore, once the experiment starts to measure sin S, then S23 can be determined with 
an uncertainty which is 6 times smaller. Thus, if a precision measurement of sin 5 is made, 
then the precision on s 2 3 becomes competitive with the precision from the disappearance 
measurement, and in fact far exceeds it at around the maximal mixing. 

Notice that by using the appearance channels we are free from the #23 octant degeneracy 
inherent in the disappearance channel. The expression of errors of s 2 3 and 5 in more general 
setting, off VOM and in matter, is given in Appendix El 

If one estimates the s| 3 error by using the simulation done by the Hyper-K collaboration 



13j, assuming 1.5 years of neutrino and 3.5 years of antineutrino runs with beam power of 



1.65 MW, one finds find that the precision on S23 exceeds the precision from the disappear- 
ance measurement around 5 — |. Of course, a full stimulation is needed to confirm this 
result, which is beyond the scope of this paper. 



C. Effect of uncertainties of #23 and # 13 on 5 

In this subsection, we discuss relative importance between uncertainties of 623 and 613 
on S. Our discussion here will be valid both on and off VOM. As already mentioned, the 
623 and the #13 intrinsic degeneracies have many common properties. For example, in near 
vacuum environment, both the #13 and the 623 intrinsic degeneracies produce 5 <H- 71 — 5 
ambiguity, and the merging of the true and clone allowed regions around \5\ = | results in 
enhanced errors of 5. Now, one can show, by the similar argument as we have given at the 
end of Sec. [TTJ that the uncertainties of 613 can affect the sensitivity to 5. Then, the natural 
question is; Which uncertainties are more important, the one of #13, or of #23? 

We argue that one can make a rough estimate of the relative importance of #13 and #23 
uncertainties on 5 determination in the following way. Notice that both of the uncertainties 
produce shift of the ellipses, up and down, along the diagonal line in the bi-probability plot 
(see Fig. [1]), and this movement results in uncertainty in 5. The uncertainties on #13 and 
#23, A#i3 and A#23, can produce the same amount of shift in the center of ellipse (the first 
term in the probabilities in (j3])) if they satisfy the relationship 

A# 

A0 23 U^ « -T-jT- (13) 
sin 6>i3 

where we have used the approximation sin2# 2 3 = 1- Using the current data of #13, 
sin 2 26*23 = 0.089 ± 0.011 Q, the uncertainty on 9 13 is 0.55°, thus using Eq. (JT3]), we can 
obtain the equivalent shift in 823, A#23| e quaz ~ 3.7°. This is comparable, or smaller than 
the current uncertainties on 623 given by 0.391 < sin 2 #23 < 0.691 at 90% CL (l4|, which 
implies A# 2 3 {current ~ 6.7° (3.8°) for upward (downward) error at lcr CL. Therefore, even in 
the current situation the uncertainty of 623 can be larger by a factor of ~ 2 than that due to 
the #13 error. Since the error of #13 is expected to become much smaller by the continuing 
#13 measurements at the time when measurement of 5 is performed, reducing the #23 error 
is of crucial importance for an accurate measurement of 5. 



IV. SIMULTANEOUS MEASUREMENT OF 6> 23 AND 5: GENERAL CASE 

In this section we present solutions of S23 and S for a given set of measurement of P and 
P at generic value of A. Instead of giving the explicit solutions we describe an iterative 
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procedure to obtain the solutions of the degeneracy equations. The explicit solutions of the 
degeneracy equations are given in Appendix 



A. Properties of the intrinsic #23 — 6 degeneracy 

To display characteristic properties of the intrinsic 823 — 6 degeneracy we derive approx- 
imate relationship between the degenerate solutions. We define differences between two 
degenerate solutions as As| 3 = (s| 3 )i — (^3)2, A(sin5) = sin^x — sin<5 2 , and A(cos5) = 
cos5i — cos<5 2 . Using 03]) we obtain 



A(cos5) cos A 



(^4® + ^e) As 23 



2 (cAq) sin 26 : 



2 

1 



2.3 



A(sin5)sinA = _ ^~^) _^fs (14) 
V 1 2 (eA Q ) sin 2^3 

under the approximation that sin 2#23 for the two solution is nearly equal, which is in fact 
an excellent approximation. This feature allows us to derive a simple relation for A (sin 6) 
as in (TH|) . which is different from the #13 intrinsic degeneracy. 

From (TH|) one can infer two important characteristic properties of the #23 — 5 degeneracy. 

• The difference in sin S comes from the matter effect, and hence A (sin 5) has to vanish 
in vacuum and the two solutions exhibit the 5 -H- tt — 5 degeneracy. 

• At any VOM, (cos A = 0), As^ has to vanish. It is because the ellipses have zero 
width at VOM and again the two solutions exhibit the 5 -H- it — 5 degeneracy. 

Therefore, a measurement near VOM can primarily determine only sin 5, and in that case 
5 -H- 71 — 5 degeneracy remains. Spectral information around VOM, if available, can be used 
to break this degeneracy. 



B. Solution to the #23 — S degeneracy equation 

We start from (JH1) and ([9]) which we call the zeroth-order solutions. Notice that we keep 
value of A in A e and A e generic here, though it has been set to | when we have discussed 
the solutions at VOM in the previous section. Under this understanding, we denote the 
zeroth-order solutions in OH]) and §§§ as (5^3)0 and 5 , respectively, in this section. By going 
back to the original equations (j3J) and using the zeroth-order solutions (]8]) and (]9]) one can 
derive the equations 

2eA 

g 23 = (^23)0 ~ sin26> 23 7- r COS A COS 5 (15) 

+ /i e ) 

(A® + Aq) sin A sin S — (A® — A ffi ) cos A cos 5 = (A® + A ffi ) sin <5 (16) 
which can be solved iteratively 4 to find the solution at generic values of A. 



4 Start with = 1/2 and calculate sin<5o using Eq. ([9]). Obtain S from the solutions of Eq. (|16|) and then 
use Eq. (|15[) to obtain s| 3 and iterate. Since sin<5o depends on sin 2#23, the change in 5 will be small and 
the iteration will converge rapidly near maximal mixing. 



s 



The second equation in (fISl) can be solved to produce the two solutions 



1 + arcsin 



(-4® + ^e) sin 5 
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where 



cos i 



(A© + A®) sin A 



V^© + A| - 2A©A© cos2A' 



• ^ - (A© - A©) cos A 

sm0 = — (18) 
v/A| + A| - 2A©A e cos 2 A 

Note, at VOM = and thus the two solutions exhibit the 5 <H- 7r — 5 degeneracy and 
they have the same value of s| 3 since cos A = 0. In vacuum, Eq. (TIE}) reduces to sin 5 = 
sin^o/ sin A, so that sin 5 is determined but cos<5 has two signs so that different values of 
s| 3 are possible, provided that we are not at VOM. 




FIG. 2: Equi-sin^ contours are plotted on ^ vs. cos A plane assuming the normal mass hierarchy. 

The typical values of ^ and cos A to be examined in Sec. |VJ are also shown by asterisks; LBNE± 
denote the cases with 20% higher (+) and lower (— ) energies than the one of the LBNE setting at 
VOM, and NuFlO and NuF15 stand for neutrino factory setting with energy 10 GeV and 15 GeV, 
respectively. 

We plot in Fig. [5] the contours of equal sin0 line on vs. — plane, assuming the normal 
mass hierarchy for which < 4^ < 1. For the inverted mass hierarchy the abscissa must be 
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interpreted as ^-. 5 In Fig. [2] we also show by asterisks the points to be examined in Sec. [V] 
with typical values of the various settings of experiments proposed. It is remarkable to see 
that in all cases sin0 is small, despite one's naive expectation that might expect become 
large when matter effects are large e.g., in the LBNE and neutrino factory settings. 

V. FEATURES OF SIMULTANEOUS MEASUREMENT OF 9 23 AND 5 IN SOME 
EXPERIMENTAL SETTINGS 

In this section, we examine features of simultaneous measurement of 623 and S by consid- 
ering several experimental settings with their typical parameters (baseline and energy). In 
the course of analysis, we illuminate some characteristic features of the 623 intrinsic degen- 
eracy. Though our analysis is based on assumed measurement of oscillation probabilities P 
and P with errors, it may capture main features of the more extensive analysis with cross 
sections, efficiencies, energy spectra etc. with systematic errors. We hope to come back to 
such an analysis in the future. We start by describing our analysis method. 



A. Analysis method 

We define the model \ 2 which is used in our analysis as follows: 

a/ fl x\ (PjM-Pjg^l (P(fl 23 ,<5)-Pfe^)) 2 
X (623,6) = — + _ (19) 

where #f 3 and S x are the test values for $23 and 8. In our analysis \ 2 is computed only at 
some fixed values of neutrino energy and spectrum information is not taken into account. 

We determine the errors to be placed on measurement of P and P in the follow- 
ing way. The number of event N in a single "energy bin" may be calculated as N = 
n T J dEf(E)^P ps rirfo^P, where rix denotes target number density and f = f(E)AE 
is the neutrino flux in a particular energy bin of width AE. Ignoring the systematic un- 
certainties on flux, cross sections, and target number density, and assuming the remaining 
uncertainties are of statistical nature only, then = oc l/y/P since AiV varies as y/~N. 
Thus, we will assume the fractional uncertainty on the oscillation probabilities are given by 



AP Pn AP 



P 





(20) 



The reference probabilities are taken as Pq = Pq = sin 2 623 sin 2 2#i 3 = 0.5 x 0.09 = 0.045, so 
that about 1100 events are assumed at the reference point, Equation ( 1201) will be referred to 
as "3% error" throughout this section. (In Sec. IV CI we will use 1% error just by replacing 
0.03 by 0.01 in (|20|) in a neutrino factory analysis.) 



In a very limited phase space can be negative. But since it does not occur in most of the proposed 
settings we do not enter into this possibility here. 
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B. Hyper-K vs. LBNE 



Let us start by examining the Hyper-K and the LBNE settings at the VOM, by which 
we mean E = 0.6 GeV and L = 295 km for Hyper-K [l3|, and E = 2.6 GeV and L = 1300 
km for the LBNE settings [l5j], respectively. In Fig. [3J the allowed contours at la, 2a, and 
3cr CL (2 d.o.f.) in 5/n vs. sin 2 9 2 3 space are presented in which a 3% error is assumed for 
both P and P measurement. Three sets of values of 5 and sin 2 #23 are used as inputs: (0°, 
0.5) shown by black contour, (45°, 0.45) by blue contours, and (—60° 0.55) by red contours. 
They are shown by + symbol in the figure, the style we keep throughout this section. 



HyperK setting: E= 0.60 GeV and L=295 km 



LBNE setting: E= 2.6 GeV and L=1300 km 
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FIG. 3: Contours of allowed regions at la, 2a, and 3a CL (2 d.o.f.) in 6/tt vs. sin 2 #23 space 
assuming 3% error for both P and P measurement. The left and right panels are for the Hyper-K 
and the LBNE settings, respectively. The true values of (S, sin 2 623) are taken as (0°, 0.5) the case 
shown in black, (45°, 0.45) in blue, and (—60° 0.55) in red. 

The features of the allowed regions shown in the left and right panels of Fig. [3] are 
remarkably similar with each other despite the larger matter effect in the LBNE setting. 
One may tempted to interpret this feature as simply due to that they are both at the VOM. 
But, actually it is not the whole story as explained shortly below. 

We observe in Fig. [3] for any one of the inputs two allowed regions corresponding to 
the true and the clone solutions, a manifestation of the 623 intrinsic degeneracy. They are 
related by symmetry under 5 -H- 71 — 5 as it is clearly visible in the figure. The symmetry 
reflects the fact that the setting is almost at VOM at which only sin 5 survives in the 
oscillation probabilities. It is also notable that true and clone contours are horizontally 
spaced, As 2 3 ~ 0, the feature which also arises by sitting at VOM, as seen in ( !T4|) . The 
allowed regions become larger as sin 2 623 get larger because the errors of P and P increases 
by rj 10% from sin 2 6» 23 = 0.45 (blue) to 0.55 (red). 

In Fig. HI the similar contours of allowed regions are plotted with the same Hyper-K and 
LBNE settings, but with higher (shown in blue) and lower (shown in red) energies by 10% 
(20%) for Hyper-K (LBNE) than their respective energies at VOM. Our choice of larger 
variations of energy off VOM reflects the fact that LBNE is equipped with a wide-band 
neutrino beam. In fact, if we choose the same rate of variation of energy for Hyper-K and 
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HyperK setting: Energy +/-10% from VOM 



LBNE setting: Energy +/-20X from VOM 
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FIG. 4: The same as in Fig. [3] but with higher (A < 5, shown in blue) or lower (A > 5, shown in 
red) neutrino energies by ±10% (Hyper-K) and ±20% (LBNE) than the ones of VOM. The true 
values of (5, sin 2 #23) are taken as (45°, 0.45) shown in blue, and (—60° 0.55) shown in red. 



LBNE the features of the contours are almost the same. The matter effect plays little role 
for this plot. It is because the slope A f°i 5 in Ifl4l is insensitive to the matter effect; it cancels 
out to leading order in A a - 2 . We observe that the direction of shift of cos 5 is opposite as 



that of s| 3 for A < | (shown in blue), while it goes along the same direction for A > | 



(shown in red), as it should be according to <^£ 

In Fig. HI we observe that while the allowed regions for the true and the clone solutions 
remain separated at la CL (apart from the red one with input = 0.55 for LBNE) 
they start to merge together at higher CL. Once the merging takes place the errors are 
governed primarily by distance between the true and the clone solutions, in addition to the 
uncertainties inherent to the measurement. If this phenomenon occurs in sizable fraction 
of the parameter space of the analysis the accuracies of measurement of both s\ z and sin 5 
are significantly lost. To our knowledge, this feature does not appear to be addressed in 
previously performed sensitivity studies for 5, while it might have been buried into the 
necessarily intricate simulation procedure. 

We notice that the movement of the allowed regions are rather dynamic despite the 
small shift of the energies in both cases. It suggests that the spectrum informations above 
and below VOM in the off-axis (narrow band) beam would be important to resolve the 
#23 intrinsic degeneracy. Naturally, the distances between the degenerate solutions become 
farther apart for larger energy variation as shown in Fig. HI suggesting that LBNE with the 
wide-band beam would be more powerful in resolving the #23 intrinsic degeneracy. 



C. A neutrino factory setting 

In Fig. [5J the similar allowed contours are presented for a particular setting for neutrino 
factory |l6] with assumed baseline of L = 3000 km. Two values of neutrino energy, E — 10 
GeV and E — 15 GeV, are chosen as typical ones which are shown in red and in blue, 
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respectively. Unlike other cases we take 1% errors for P and P, which implies ~ 10 4 events 
in each channel at the reference probabilities Pq = Pq = 0.045. Despite the smaller errors 
taken discriminating power of the clone solution from the true one is quite limited. Two 
regions of true and clone solutions merge apart from those in region of sin 5 < at E = 10 
GeV. Since they are below VOM (A < |) the relationship between true and clone solutions 
are similar at both energies. 

NuFactory: L=3000 km, E=10, 15 GeV 

0.65 
0.60 
0.55 

CO 

c\! 

1> 

CM 

B 0.50 
0.45 
0.40 
0.35 

-1.0 -0.5 0.0 0.5 1.0 

6/ir 

FIG. 5: Contours of allowed regions at la, 2a, and 3a CL (2 d.o.f.) in S/tt vs. sin 2 #23 space for 
neutrino factory setting with L = 3000 km. Two particular values of neutrino energies are taken 
for comparison, one at E = 10 GeV (shown in red) and E = 15 GeV (shown in blue). Unlike other 
cases 1% error is assumed for both P and P measurement. The true values of (5, sin 2 #23) are 
taken as (45°, 0.45) and (-60° 0.55). 




VI. A SETTING AT THE SECOND OSCILLATION MAXIMUM 

In this section we discuss some new aspects that show up when the beam energy is tuned 
to, or it covers over, the second VOM, A = The most important change in this region 
is that A & cx A becomes larger by a factor of three. Then, S dependent term which comes 
from interference between A§ and A & gets larger by the same amount. It influences to the 
errors of s| 3 and 5. As 2 ^ in ( fTTj) remains unchanged if the uncertainties AP and AP are the 
same at the 2nd VOM. On the other hand, A(sin<5) becomes smaller by a factor of three, 
implying a better sensitivity to CP at the 2nd VOM. Therefore, Eq. (Tl2|) is modified to 
A(4 3 )~|A(sin5). 

To elucidate the new features of the setting at 2nd VOM we examine the case with ESS 
(European Spallation Source) [l7| , the facility to be build in Sweden which plans to produce 
5 MW of beam power which can be used to construct intense neutrino beam. The baseline 
to a detector is assumed to be L = 540 km and the energy at around the 2nd VOM is given 
approximately as E = 0.353 GeV. In Fig. O presented are the contours of allowed regions by 
assuming 3% measurement of P and P. We take neutrino energies 3% higher (0.364 GeV, 
A < |, shown in blue), or 3% lower (0.342 GeV, A > |, shown in red) than that of 2nd 
VOM. 
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ESS setting: Varying E by +/- 3% 
0.65 ' 
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FIG. 6: Contours of allowed regions at la, 2a, and 3a CL (2 d.o.f.) in 5 /it vs. sin 2 623 space for 
ESS setting with L = 540 km. Two particular values of neutrino energies are taken for comparison, 
one at E = 0.342 GeV (shown in red) and E = 0.364 GeV (shown in blue), which amount to ±3% 
of the one at 2nd VOM. The true values of (5, sin 2 ^) are taken as (45°, 0.45) and (-60° 0.55). 

One notices that the allowed regions significantly shrink along S direction, the effect we 
mentioned above. It is also remarkable that clone solutions moves significantly upward or 
downward despite very small variation of neutrino energy of ± 3%. It implies that resolving 
power of the 823 intrinsic degeneracy, in principle, can be higher at the 2nd VOM. 

VII. CONCLUSION 

In this paper, we have shown that S23 can be determined automatically at higher accura- 
cies than sin S in experiments which can measure S using v e and v e appearance channels. The 
accuracy for s| 3 determination supersedes that of the conventional method which relies on 

disappearance channels in region 40° < 823 < 50°, assuming high enough statistics that 
allows reasonable precision on the measurement of CP phase, S. The reasons for such high 
accuracy is that the appearance method is free from the problems of Jacobian broadening 
in the translation sin 2 2#23 — > for values of 623 near 45°, and of the determination of the 
octant of 6*23- 

With use of v e and v e appearance channels, determination of Q23 is inherently coupled 
with measurement of CP phase 5. Thus, the uncertainties on the measurement of 623 and 
S are correlated. Furthermore, our treatment reveals that measurement of the two indepen- 
dent parameters metamorphoses into a new strategy which may be called as "simultaneous 
measurement of 623 and <5" . Reflecting on the nature of these measurements, a new parame- 
ter degeneracy called the "6*23 intrinsic degeneracy" is identified. In general, two degenerate 
solutions of #23 and are given for each set of measurements of oscillation probabilities, P 
and P, producing a new ambiguity and a new source of uncertainties, hitherto unrecognized. 
This degeneracy must be resolved to make an accurate determination of both 623 and S. 

To understand characteristic features of correlated measurement of 823 and 5, we have 
undertaken the following two approaches: (1) Toward resolution of the degeneracy, we have 
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analyzed the natures of the 823 intrinsic degeneracy, and obtained the analytic formulas for 
the degeneracy solutions which help us to understand characteristic features of the degen- 
eracy. (2) We have simulated simple data set of P and P with suitably assumed errors 
to obtain the contours of allowed regions in S23 vs. 5 space. We have examined various 
settings, Hyper-K, LBNE, and neutrino factory all around (or below) 1st VOM, and ESS 
(European Spallation Source) near 2nd VOM, using their typical values of baseline and 
neutrino energies. 

From these studies we have observed the following features of our new method, simul- 
taneous measurement of 6*23 and 5. The obtained contours of allowed regions confirm our 
expectation that the error of S23 determination is approximately six times smaller than that 
of sin 5 apart from the problem of merging the degeneracy solutions. Importantly, the move- 
ment of clone solutions against varying energy is generally quite large, and the effect is more 
significant around 2nd VOM than the 1st VOM. Thus, we have uncovered the crucial point 
that spectrum information is of key importance to resolve the 823 intrinsic degeneracy. This 
feature is reminiscent of the one encountered in the #13 intrinsic degeneracy. 

In addition, the contours of allowed regions revealed some characteristic features of the #23 
intrinsic degeneracy, such as 5 — > ir — 5 symmetry, in most of the settings we discussed. While 
the feature must hold at VOM or in vacuum, its robustness in fact requires explanation. 
We have shown that it can be understood by our analytical expressions of the degeneracy 
solutions, which implies that the matter effect plays only a minor role in the relationship 
between the two solutions of S. 

We have concluded that our new method, using u e appearance channels, offers a better 
way of measuring 823 with high accuracy, in the region 40° < 823 < 50° in particular, and 
the simultaneous measurement of 823 and 5 is the right strategy to carry out precision 
measurement of the both parameters. Incorporation of these ideas into a full stimulation 
will be undertaken in a follow up paper (l8| . 
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Appendix A: Uncertainty of 5 from Probabilities and s\ 



In this Appendix, we calculate the uncertainties in cos 5 and sin 5 which would be pro- 
duced by small measurement uncertainties in P, P, and s 2 ,^- Here, we do not consider as 
a parameter determined by our appearance method with (P, P) measurement, but regard it 
as an external parameter determined e.g., by using the disappearance channel. 

If we start with the appearance probability expressions, Eq. 01]), and expand the cos(A±<5) 
terms: there are two terms with 5 dependence, one proportional to sin 5 and the other 
proportional to cos 5. We then use the chain rule of differentiation, to leading order in 

^23 <^ 1, to obtain 6 



e = 



and for cos 2# 2 



A (cos 5 cos A) 
A(sin<5sin A) 



tan# 
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eA Q 
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eA Q 



AP 



+ 



A ffi ± An 



eA P 
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'23 



(Al) 



We interpret this result as follows: the uncertainties on the measurement of (sin S sin A) or 
of (cos S cos A), assuming that the uncertainties on P, P, and are uncorrelated, are given 
by 



A (cos S cos A) 
A (sin 5 sin A) 




(VP) 2 + 



A, 



eA 







(vpy 



A m ± Ac 



eA 







(Vs 2 ^ 



'23^ 



(A2) 



where VO = ^ denotes fractional error of the quantity O. 

In vacuum, when A® = A®, it's clear that uncertainty on sin 5 does not depend on the 
uncertainty of whereas the uncertainty on cos 5 depends strongly on this uncertainty. 
In matter, however, the uncertainty on both sin 5 and cos 5 will depend on the uncertainty 
on S23. One needs both measurements to make a precise measurement of 5 otherwise the 
measurement will suffer from a degeneracy similar to the 5 <H- 7r — 5 degeneracy seen at 
VOM. In the Introduction we simplified this rather complicate situation to make a point 
regarding the importance of an accurate determination of S23. 



In near vacuum environment, 



Am 31 



< 1, Afi w A ffi and we have t -^piA 



eAr. 



1.5, 



which entails the coefficients in y/^ x (5X) 2 in ( 1A2j) of order unity. Therefore, in near 
vacuum, an estimation of the uncertainties on (sin S sin A) (or (cos S cos A)) are given by a 
very simple expressions 



Uncertainty of (sin S sin A) ~ f c 
Uncertainty of (cos 5 cos A) ~ f c 



(VP) 2 + (VP) 2 , 



(VP) 2 + (VP) 2 + 4(Vs 2 3 ) 2 



(A3) 
(A4) 



where /error is a coefficient of order unity. 7 Thus, for | sin 51 near 1, where a precise measure- 
ment of cos S is important for an accurate determination of the angle S, the uncertainty in 



6 Eqs. (|XTJ) and $E2j) is more accurate if you interpret P and P in the denominator in the hrst two terms 
of these equations as (P) and (P), the ones averaged over S. In other words (P) and (P) imply the 
coordinate of the center of the bi-P ellipses. 

7 In matter there exists 77(1? s| 3 ) 2 term in the square root in (|A3I) but with small coefficient. For example, 
T) ~ 0.24 for LBNE. 
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sf 3 plays a significant role. The feature that sin 5 does not depend on s| 3 while cos S does in 
the near vacuum environment can be easily understood by the bi-probability plot given in 
the right panel of Fig. HJ When S23 is varied within the uncertainty "center of gravity" of 
the ellipse moves along the diagonal line in P — P plane, which affects predominantly cos 8 
but not so much to sin 5, as this is a measure of the distance away from the diagonal. 

The results in (1A3I) and (1A4I) are consistent with each other because here sin 5 and cos 5 
are treated as independent. Otherwise, the relationship A(cos<5) = — tan<5A(sin5) must 
hold. 



Appendix B: Errors of S23 an d 8 °ff the VOM 

One may ask the question; Can one derive the formula of errors of S23 and 5 analogous 
to pip in general setting off VOM and without vacuum dominated approximation? In this 
Appendix we give an answer to this question. 

We start from the expressions of P and P in (j3J). We denote the changes of P and P 
induced by a small changes of S23 and 5 as AP and AP. Conversely, it may interpreted as 
the relations between small shifts of S23 and S and their cause, i.e., errors of P and P. To 
leading order in the small variations we obtain the relationship between the shifts As| 3 , AS, 
AP and AP as 



AP 4^ si n (A + 5) + AP 



sin(A - 6) 



2 v^4e^e [ (^e + -4®) sin A cos S + (A© - A©) cos A sin 6] 



(Bl) 



A5 



ap (£) - ap (I; 



2sin26 l 2 3 (eA ) [(A® + A©) sin A cos S + (A© - A©) cos A sin 5] 



(B2) 



where we have ignored the variation of sin 2# 2 3 to derive these expressions, since we are only 
interested in the region near maximal mixing. 

From these expressions it is straightforward to derive the formulas of errors of S23 and S 
by assuming that the errors of P and P are independent with each other, and by ignoring 
errors in neutrino energies, Am| 2) and etc. They generalize ffTTl) to the case in matter and 
off VOM. Notice that these expressions of errors are obtained in our appearance method for 
simultaneous measurement of #23 and 5, and they should not be confused with the errors 
discussed in Appendix [S] 



Appendix C: Solution of the Degeneracy Equation 

Here, we obtain explicit solutions to the degeneracy problem involving S23 and S. 
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1. General solutions within the normal hierarchy 

By using (j3J) we obtain 

cos A cos 5 
sin A sin 5 
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sin 2#23^4© 
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(CI) 



Then, the equality cos 5 + sin 5 = 1 can be written in a form of quadratic equation for 



x = s 



2 . 
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(C2) 



where 



(A% + A%) - 2A ffi /4 ffi cos 2A + Ae 2 A 2 Q sin 2 2A 

+ -^A e ) cos 2 A + At 2 A% 
A ® A © / 
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3 + (S - 2 (£)(£ lcos2A 



(C3) 



The solution to (1C2|) is of course given by 



x 



b±Vb 2 - 4ac 
2a 



(C4) 



The solutions for cos 5 and sin 5 can be obtained by inserting (1C4j) into (IClj) . The discrimi- 
nant D = b 2 — 4ac is given by 



P = sin 2 2A 



+ 4 £ Mj{ 2 (P + P>- £\£) - 2c08 2A(£, ffi + £ie-££)} 



+ 16e 4 ^ sin 2 2A 



(C5) 



When D > (1C4I) describes two solutions of S23 and 5 for a given set of observable P and 
P. By inserting them into fICll) one obtains the corresponding solutions of cos 5 and sin 5. 

The degeneracy solutions within the alternative (flipped A-sign) mass hierarchy can be 
obtained by making replacements A — > — A, A 9 — > —A e , and A & — > — A^ in flC4j) and (IC3I) . 



2. Perturbation around the zeroth-order solution 

We obtain perturbative solution of S23 = a; and 8 in a form x = x + x%, and 8 = 8 + 5\ 
to leading order of small deviation of A = | — y. Here, xq and 80 imply the zeroth-order 
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solution given in flSJ) and Q). 8 Notice that cos A m y, cos2A rs — 1 + y, and sin2A 2y. 
We obtain the equation for x\\ 

(A® + A e ) x 2 

= 4e 2 ^ sin 29 23 y 2 + 4A e A e x^ 2 - 2 ( £-A® + ) x o2 / 2 + ^~^~y 2 (C6) 

where a small changes to sin 2#23 due to correction to Xq is ignored. The solution to this 
equation is given by 

y 

Xx = ±jj — ^-2eA sin 2^23 cos <5 (C7) 

where we have used Using sin(5 + Si) ~ sin<5 + Si cosS in (ICip one obtains, ignoring 
order y 2 term, 

A e — A® A® — A® 

2eA sm 26*23 ^4® + ^4® 



where we have used (jC7j) . The correction to 5 comes from the matter effect, whereas the 

2 

23 



one to dominantly vacuum effect 
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